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Abstract 

Fullerenes have several advantages as potential materials for organic spintronics. Through a 
theoretical first-principles stndy, we report that fullerene Ceo adsorption can induce a magnetic 
reconstruction in a Ni(lll) surface and expose the merits of the reconstructed C6o/Ni(lll) spinter- 
face for molecular spintronics applications. Snrface reconstrnction drastically modifies the magnetic 
properties at both sides of the Cgo/Ni interface. Three outstanding properties of the reconstructed 
structure are revealed, which originate from reconstrnction enhanced spin-split tt— d coupling be¬ 
tween Cqo and Ni(lll): 1) the Cgo spin polarization and conductance around the Fermi level are 
enhanced simultaneously, which can be important for read-head sensor miniaturization; 2) local¬ 
ized spin-polarized states appear in Cgo with a spin-filter functionality, and 3) magnetocrystalline 
anisotropic energy and exchange coupling in the outermost Ni layer are reduced enormously. Sur¬ 
face reconstruction can be realized simply by controlling the annealing temperature in experiments. 
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Non-magnetic/ferromagnetic interfaces of magnetic hard disk drives are crucial for their 
performance. To achieve read-head sensor miniaturization, it is essential to minimize the 
resistance of the system to maintain an ideal data-transfer rate and signal-to-noise ratio; 
at the same time, a large magnetoresistance (MR) ratio is decisive for optimum function¬ 
ality of read-heads[l, 2], It is therefore desired to hnd a method that can simultaneously 
increase the MR ratio and the conductance of the device j^. Recently, tuning the properties 
at organic/ferromagnetic interfaces by aromatic molecules has attracted broad attentionP- 
Due to the tt — d hybridization, chemical adsorption of aromatic molecules on magneRc 
surfaces produces new spin-split hybridized states at the interface (called a spinterface 


Q. 






These states can be used to produce thermal-robust molecular spintronic devices [l7(]. A 
challenge of this method is that chemical adsorption usually broadens the molecular orbitals 
near the metal Fermi level, which acts against the desired appearance of spin-polarized and 
energy-concentrated states [l7, 181. Another possible usage of the organic/ferromagnetic 
interface is that the effective magnetic coupling of surface atoms can be modihed by the 


adsorbate 


3, 19 


2l| . Thus it is possible to create hard/soft composite magnetic struc¬ 


tures by self-assembly in molecular adsorption to achieve desired applications in permanent 
magnets, recording media and spintronics 22 ]. 


Fullerenes anc 


ganic devices [23 


their derivatives are building blocks of potential high performance or- 


281]. Meanwhile, it has been proven that Cgo adsorption can induce non¬ 


magnetic metal surface reconstruction, i.e. rearrangement with different bonding of surface 


atoms 


29, 


These reconstructions have decisive influences on the their charge transport 


properties [Mi- Thus, extending to magnetic metal surfaces, one can expect that adsorption- 
induced reconstruction could also have signihcant effects on the spin transport properties 
of fullerene/ferromagnetic interface. It is crucial to identify the existence and magnetic ef¬ 
fects of reconstruction at these organic/ferromagnetic interfaces. For C6o/Ni(lll), height 
prohle measurements with scanning tunnelin g ni icroscopy (STM) show that Cgo can adsorb 
at different heights above the Ni(lll) surface 32j, which is a hint for surface reconstruction. 
However, the reconstructed atomic structure and the effects of reconstruction on interface 
magnetic properties has not been studied to our knowledge. 

In this letter we investigate the geometric and magnetic properties of the C6o/Ni(lll) in- 
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terface by first-principles methods. We show that the reconstructed structure is energetically 
favored over the unreconstructed one. We demonstrate that, in comparison with the unre¬ 
constructed structure, the reconstructed one has the following superior properties: 1) the 
density of states (DOS) and the spin polarization of Cgo are enhanced simultaneously around 
the Fermi level; 2) the molecular spin-polarized states are concentrated in energy around the 
Fermi level; 3) the magnetic coupling and magnetocrystalline anisotropic energies (MAE) 
of atoms in the outermost substrate layer are signihcantly reduced. The above changes in 
properties show that one can signihcantly affect magnetism at the organic/ferromagnetic 
interface through surface reconstruction. This prediction could have further applications in 
molecular and organic spintronics. 

Calculations were performed using the plane-wave-basis-set Vienna ab-initio simulation 
package (VASP) 33(]. The Ni(lll) surface was modeled by a hve-layer-slab with a 4x4 sur 
unit cell per Cgo, which cell size is determined from electron diffraction experiments 34j . 


ace 


Projector augmented wave potentials 


35|. 


were employed with a kinetic energy cutoff of 500 


eV and with a iP-point sampling of 4x4. For the exchange-correlation functional, the 
Perdew-Burke-Ernzerh of generalized gradient approximation was utilized 37|. The calcu¬ 
lations of magnetic couplings and MAE utilized the Quantum Espresso package j^| with 
equivalent computational parameters as in the VASP calculations. The adsorption struc¬ 
tures with and without reconstruction were selected by (a) symmetry- and size-matching 
between Csn and the Ni(lll) surface, and (b) the informations of the same and similar 


systems 


32, 39 


43| (details see the Supplemental Material j44|j. 


Fig. [T] shows the most stable structures for both the reconstructed and unreconstructed 
cases (Rec and Unrec, resp.). Their stabilities, i.e., adsorption energies corrected for va¬ 
cancy formation energies, are —2.13 eV for Unrec and —3.41 eV for Rec (see Supplemental 


Material 


M). 


As shown in Fig. [T^ and b, Cgo binds with Ni through a Cg hexagon parallel to the Ni(lll) 
surface. In Unrec, the centre of gravity of Cgo locates at a bridge site of the surface. Six 
C—Ni bonds are formed with bond lengths ranging from 2.0 to 2.05 A. In Rec with a 7-atom- 
cavity (Fig. [T]d), the molecule locates at an fcc-hollow site of the second Ni layer. Three C 
atoms of the Cg hexagon bind to three Ni( ”q” in Fig. [TJl) in the second layer and six C in 
the second layer of Cgg bind to six Ni(”o”) atoms in the top surface layer, forming nine C—Ni 
bonds in all. The corresponding bond lengths rang from 1.93 to 2.00 A, indicating stronger 
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FIG. 1. (a) and (b): Top views of the unreconstructed (Unrec) and reconstructed (Rec) structures 
of C6o/Ni(lll); only the bottom part of Cgo and the top layer(s) of Ni are shown; the dashed lines 
outline a Ni(lll)-(4x4) surface unit cell. Different types of surface atoms are labeled with different 
colors classified by their magnetic moments and their distances from the Cgo carbon atoms, (c) 
and (d): Magnetic moments (in /tb) of surface atoms in Unrec and Rec structures. 


bonding strength than the one in Unrec. The formation of more and stronger C—Ni bonds 
in Rec snrpasses the energy cost of forming a seven-atom hole in the snrface and makes the 
Rec structure more stable. The C6o/Ni(lll) system was investigated previously by STM 


32|, and two adsorption configurations were fonnd. Ceo in one configuration was fonnd to be 


2.2 A lower than the other from the measured apparent height profile 321. The Ceo height 
difference between Rec and Unrec in onr calculation is 2.0 A, by defining the molecular 
height as the height of the outermost hexagon, and using the average height of the second 
snbstrate layer’s Ni atoms as a reference. Thus, the low and high configurations in the STM 
experiment can be interpreted to be the Rec and Unrec structures in onr calculation. The 
coincidence in molecular height between theory and experiment supports the reliability of 
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our proposed structures. 

We quantify the spin-dependent charge transfer by the Bader charge analysis method 


45|. 


In Unrec, 0.64 electrons are donated into each Cgo from the metal substrate, among which 
0.26 electrons occupy spin up states while the other 0.38 occupy spin down states. The 
charge transfer increases to 1.82 electrons in Rec, with 0.85 electrons occupying spin up and 
0.97 electrons occupying spin down. These charge transfer have significant effects on the 
electronic and magnetic properties of both the adsorbed Cgo and the Ni surface, as detailed 
below. 

We plot the spin-polarized projected density of states (PDOS) of Ceo for Rec and Unrec 
in Fig. [2^, and the corresponding spin polarization ratio (SPR) in Fig. [2 ]d. The SPR is 
defined as SPR(U) = (DOSt(^) - DOS;(U))/(DOSt(^) + DOS;(U)). Fig. EJr shows that, 
compared to the free Cgo monolayer, the interaction with Ni broadens the molecular orbitals 
due to the hybridization of Cgo P orbitals and substrate d states llj] . These hybridization 
enhance the PDOS of Ceo near the Fermi level, changing the molecule from semiconducting to 
metallic. In particular, the PDOS of Cgo at the Fermi level in Rec is 2.5 times that in Unrec. 
Therefore, we can expect a higher conductance in Rec than in Unrec. More importantly, 
due to the magnetic surface, the molecular DOS is spin polarized after adsorption, especially 
near Fermi energy. According to the Julliere model of spin-dependent tunneling^, the 
MR of the system is positively correlated with the SPR at the Fermi energy. From Fig. 
|2 ]d, the SPR in Rec is about 19% around the Fermi energy, almost three times that in 
Unrec, which is about 7%. Therefore, from the PDOS and SPR, we can expect that the 
surface reconstruction enhances the conductance and MR of Cgo simultaneously. This feather 
meets the requirement of the miniaturization of read-head sensors, as we mentioned at the 
beginning js I. 

Another notable character of Rec is the DOS distribution in Cgo- We plot the average 
PDOS in the bottom, the middle, and the top parts of Cgo (Fig- Et)- It can be seen that 
these PDOSs are quite different from each other. These differences mean the delocalized 
molecular orbitals near the Fermi level in Cgo are broken into groups of localized orbitals. 
We plot the local DOS in the energy interval of [0, 0.2] eV in Fig. [2]i to show the orbital 
distribution (the shape of the local DOS in [—0.2,0] eV is similar). We can see that the 
states are localized at the equator of the Cgo- In particular, from the middle PDOS in Fig. 

PDOS at the Cgo equator has a rather strong intensity near the Fermi level, and that the 
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FIG. 2. (a) PDOS of free Ceo monolayer, and of Ceo in Rec and Unrec structures, (b) Spin- 
polarization ratio (SPR) of Cgo in Rec and Unrec. (c) PDOSs of different parts of Cgo in Rec: 
the bottom (blue), the middle (red), and the top (black); the colors of PDOS curves match the 
colors of the atoms in the inserted CgOj the metal surface being horizontal and below in this view; 
the green dashed lines enclose the most suitable energy range for the molecule to filter spin, (d) 
Spatial distribution of spin down local DOS in the energy interval of [0, 0.2] eV (isovalue 0.0015 
e/Bohr^), viewed along the surface. 


corresponding spin polarization at the Fermi level is significantly stronger than in the other 
two parts of the molecule. The PDOS at the equator has two advantageous features. One is 
that the Fermi level passes through the spin polarized peak, so these states can be accessed 
by a low bias voltage to reduce power dissipation. The other is that the polarized states are 
concentrated around the Fermi level so that a large current can be expected. Thus, the Rec 
system can hopefully be used as a basis of building high efficient molecular spin filters that 
work in the bias window of ±0.2 eV, as is enclosed by green dashed lines in Fig. This 
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feature is not present in the Unrec structure 28| . The reason is as follows. In Unrec, only 
0.6 electrons transfers into the three-fold degenerate lowest unoccupied molecular orbitals 
(LUMO) of Cgo- Such a small electron transfer cannot signihcantly shift the Cgo LUMO 
toward the Fermi level (Fig. |2^). On the contrary, because of the reconstruction, 1.82 
electrons transfer to Ceo, and one LUMO orbital is supposed to be occupied. The partially 
occupied LUMO is the source of the states at the equator of Cgo- The layer-dependent DOS 
character of Ceo on Ni(lll) is similar to the DOS character of a double-decker molecule after 
adsorption}]^. The difference is that here the reconstruction plays a key role. 

In addition to the above changes in the molecule, the reconstruction also signihcantly 
modihes the magnetic structures of the surface. The calculated magnetic moments of the 
outermost nickel atoms in the Unrec or Rec cases are presented in Fig. [TJ:; and [T]i. Based 
on the magnetic moment values and the internuclear distances from the molecule, the nickel 
atoms at the interface can be grouped into four types for both Rec and Unrec. For the 
unreconstructed case, four Ni atoms are strongly affected ( ”i” and ”j” in Fig. [T}::): these 
four, which are right beneath the molecule in the outermost Ni layer, hnd their magnetic 
moments reduced to about 0.3 pb- Next to these four atoms, four other Ni (”k”) atoms are 
slightly inhuenced, converting their magnetic moments to 0.63 pe- Others (i.e. ”1”) keep 
their clean surface value of 0.69 /ib. By contrast, for the reconstructed case, all magnetic 
moments of the hrst Ni layer are signihcantly reduced. Six nickels (”o”) along the rim of the 
hole change their magnetic moments to 0.28 /tb while three others (”p”) change to about 
0.44 pb- As the reconstruction removes seven top layer atoms, the Ceo directly bonds to 
three nickels of the second layer (” q”) and consequently reduces their magnetic moments to 
0.35 pb- 

To further investigate the magnetic properties of the surface, we calculated the magnetic 
coupling and MAE of the surface atoms. We used a Heisenberg model with a Hamilto¬ 
nian H = fo describe the hrst nearest neighbor magnetic interaction 

between Ni atoms and dehned the magnetic coupling strength between atom a and atom 
(3 as was Calculated from energy diherences between properly selected 

couples of magnetic conhgurations For the clean unreconstructed surface, the coupling 
is 8.5 meV. From Table H] we can see that the adsorption softens all the magnetic couplings 
between the surface atoms. This is similar with a previous study on Co interfaces 2^. The 
decreases of these magnetic coupling strengths are strongly correlated with the distance from 
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TABLE I. Magnetic coupling strengths Jai 3 fJ-afJ>i 3 (in rneV) between different types of mutually 
bonded atoms at the surface. Here a,l3 = for Unrec and a, (3 = o,p,q,r for Rec indicate 

different types of atoms as labelled in Fig. [Tt and d. 


Types(Unrec) 


Types(Rec) 


ii 

1.14 

oo 

2.16 

ij 

1.21 

op 

1.92 

ik 

2.51 

or 

2.22 

il 

2.39 

pr 

1.08 

jk 

2.62 

qq 

1.95 

jl 

2.94 

qr 

4.68 

kl 

7.21 

rr 

6.97 

11 

7.00 

- 

- 


the molecule and the bonding condition to the molecule. 

To analyze the reason of this softening, we examined the PDOS of the d orbitals of 
selected Ni atomsjd^. These orbitals are grouped into d-,,, which have out-of-surface-plane 
components {d-,, = dxz + dyz + dz^), and d^ lie within the surface plane. We find a notable 
change is that the spin down orbitals move to lower energy and the spin up orbitals move 
oppositely so that the spin-split energy {Ess) is reduced. We list the d-band-center shifts 
of the selected atoms relative to the clean unreconstructed surface and the corresponding 
decreases in Ess in Table [TTl We can see that the decrease of Eg'’^ + Eg^'s on these atoms are 
qualitatively consistent with the reduction of the magnetic moments of the corresponding 
atoms, i.e. on the order oi o ^ i > q > p. Another feature is that the spin down orbitals shift 
relatively more in the Ni atoms which have C—Ni bonds. Thus, the C—Ni bonds play an 
important role in the orbital shifts. As is demonstrated in valence bond theory Sd metals 
usually move some d electrons of majority spin into its minority orbitals and use the empty 
d orbitals to form new hybridized orbitals when forming bonds with organic compounds. So 
we conclude that the spin down orbital shifts are caused by hybridization. Meanwhile, as 
the substrate loses electrons, the center of the total d orbitals must move to higher energy. 
Therefore, the spin up orbitals move toward the Fermi level. These are the reasons for the 
relative energy shifts of d orbitals. 
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TABLE II. Energy shifts of d band centers relative to the clean unreconstructed surface and of 
Ess for Nio, Nip, Nig and Nij (in eV) defined in Eig. [TJ The reference d band centers are —1.49, 
—0.81, —1.39, and —0.74 eV for d“^, d“^, and respectively. 


Atom types 

d7 

^down 

Tpdcr 

^SS 

d7 

^down 

'crd'K 

^SS 

Nio 

0.22 

- 0.19 

0.27 

0.14 

- 0.21 

0.30 

Nip 

0.26 

- 0.01 

0.41 

0.21 

- 0.03 

0.41 

Nig 

0.08 

- 0.20 

0.40 

- 0.03 

- 0.36 

0.32 

Nh 

0.09 

- 0.28 

0.31 

- 0.01 

- 0.41 

0.25 

Clean surface 

0 

0 

0.68 

0 

0 

0.65 

The MAE are calculated from the difference of spin-orbit interaction (SOI) 

energies 

for spins along different 

axes. The spin-orbit interaction can be calculated by non- 

collinear density fun 

LCti 

one 

d theory (DFT) 

as well as a 

second order perturbation Esoi = 

\2 |(o|L-S|ii)|2 

Eu — Eo 

49 

51 

, where o and 

u denote occupied and 

unoccupied 

collinear 


Kohn—Sham orbitals and A is a coupling constan 


tained by linking the results of these two methods 4^. Values along the axes [111] (out of 


The MAE of each layer can be ob- 


the surface plane), [101] and [112] (in plane, see Fig. [H::) were obtained. The clean unre¬ 
constructed surface has an in-plane magnetization with an MAE of 0.24 meV/atom with 
degenerate energy for the easy axes [112] and [101]. After Ceo adsorption, the MAE of the 
first layer becomes 0.14 meV/atom and 0.03 meV/atom in Unrec and Rec, respectively. This 
behavior coincides with the knowledge of the origin of MAE that the enhancement of vertical 
interaction (by molecule here) can weaken the stability of in-plane magnetization 521. 

From the above results, we can see that the adsorption of Cgo on Ni(lll) will soften 
the magnetism of the outermost Ni atoms. Such softening is significantly enhanced by 
reconstruction. Therefore the reconstruction can be used as a convenient way to generate 
hard/soft composite magnetic structures to realize specihc functions 22]. 

In summary, we determined the atomic structures at the C6o/Ni(lll) with and without re¬ 
construction. The reconstruction not only stabilizes the molecule but also causes significant 
changes in the spintronic and magnetic properties at both sides of the organic/ferromagnetic 
interface. On the molecule side, by reconstruction , the DOS and SPR of the adsorbed 
molecule at Fermi level are increased, which could simultaneously improve the conductance 
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and MR of the system. The reconstruction also creates a new spin-polarized and energy¬ 
concentrated state at the equator of Ceo near the Fermi energy, which makes it possible 
for the molecule to be used as a spin hlter. This suggests that, instead of using double- 
decker molecules, one can also use three-dimensional molecules, such as fullerenes, to obtain 
layer-dependent spin-polarized states. On the surface side, reconstruction drastically re¬ 
duces the exchange coupling and MAE of the outermost layer. The change of the exchange 
coupling can be related to the d-bands shift under the influence of molecular adsorption. 
These findings reveal the importance of reconstruction on the organic/ferromagnetic inter¬ 
faces, and could serve as basis for developing novel spintronic devices. Besides, another our 
investigation found that Cgo could induce a different type of reconstruction on the ferro¬ 
magnetic Fe(lOO) 53|. Combining this with the knowledge that Cgo induces non-magnetic 
metal surfaces reconstructions in different types, it is reasonable to believe that various re¬ 
constructions can also happen at other interfaces between different magnetic surfaces and 
molecules (e.g. C 70 , thiolates, and graphene). Thus our discoveries can be extended to other 
systems with various combinations between organic materials and magnetic surfaces. 
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